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This paper presents some preliminary evaluation results of a full turbofan engine ice 
crystal icing altitude scaling application for a heavily instrumented ALF502R-5 engine test 
conducted in the NASA Glenn Research Center Propulsion Systems Laboratory (PSL) dur- 
ing October and November of 2015. This engine, S/N LF11, was instrumented with four in- 
ternal video cameras to observe various ice buildup processes downstream of the second exit 
guide vane trailing edge (i.e. EGV2 TE) region of the low pressure compressor (LPC). To 
determine possible scale test conditions at lower altitudes a reference condition chosen from 
one of this engine’s key rollback test points, i.e. FLT 850, was selected. The scaling method 
considered is based on the thermodynamic model developed previously to describe the po- 
tential ice accretion on unheated surfaces inside a turbofan engine compression system from 
fully glaciated ice crystal clouds. The PSL engine icing test capability has been demonstrat- 
ed in its 2013 inaugural ice crystal icing test. But the rather limited ice crystal cloud calibra- 
tion range available at the time and the lack of reliable measurement instrumentation need- 
ed along the internal core flow passage of the LPC has prevented one to setup the desired lo- 
cal airflow and ingested ice crystal cloud thermal and phase conditions according to the scal- 
ing method. Consequently a simplified altitude scaling approach was investigated. The 
NASA in-house engine icing risk assessment code COMDES was utilized to determine possi- 
ble scale test conditions at lower altitude with similar engine icing risk potential by best 
matching the following three icing related parameters of the reference condition: (1) the lo- 
cal air static wet bulb temperature and (2) the local ice crystal melt ratio of a prescribed 
smaller particle size at the EGV2 trailing edge location, and (3) the engine fan face total wa- 
ter content to air density ratio of the ice crystal cloud. The instantaneous images and the 
time-lapsed movies of ice buildup in the EGV2 TE region were acquired in LF11 and are 
used with other relevant air and ice crystal thermodynamic and engine performance data to 
assess the lower altitude icing simulation results. It is found from the preliminary evaluation 
that altitude scaling is possible in PSL to simulate the important ice crystal accretion fea- 
tures and the associated icing effects on engine performance. 


Nomenclature 


= Modified relative heat factor, dimensionless 

= Specific ratio of ice particle kinetic heating to latent heat absorbed from melting, dimensionless 
= Specific heat of air, cal/g K 

= Specific heat of water at the surface temperature, cal/g K 
= Convective heat-transfer coefficient, cal/s m?K 

= Gas-phase mass-transfer coefficient, g/s m? 

= Gas-phase mass-transfer coefficient, g/s m? 

= Ice water content, g/m? 

= Non-dimensional ice water content, = [WC/pa 

= Liquid water content of melted ice, g/ m> 

= Ice crystal melt ratio, = LWC,/TWC,, dimensionless 

= Local Mach number, dimensionless 

= Ice crystal median mass diameter, um 

= Water droplet median volume diameter, um 

Fan speed, rpm 

= Core speed, rpm 


* Principal Research Scientist, 21000 Brookpark Rd., MS 11-2, AIAA Associate Fellow. 
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m = mass flux of water evaporated per unit time, Ibm/ft? s 


e 


Miny » = mass flux of ice/liquid-water particle impinged per unit time, lbm/ft” s 
‘nelt.imp = ‘Mass flux of water-melt impinged per unit time, lbm/ft* s 

mo = Stagnation melting fraction, dimensionless 

no = Stagnation freezing fraction, dimensionless 

D = Pressure, N/m? 

Pw = Vapor pressure of water in atmosphere, N/m? 

Pww = Vapor pressure of water at the icing surface, N/m? 

conv = Surface heat loss due to convection, Btu/hr ft? 

evap = Surface heat loss from evaporation, Btu/hr ft? 

freeze = Surface heat gain from release of latent heat of fusion, Btu/hr ft? 

kinetic = Surface heat gain from kinetic energy of ice crystals and water drops, Btu/hr ft? 

melt = Surface heat loss from release of latent heat of fusion, Btu/hr ft? 

RH = Relative humidity, dimensionless 

RTD1 = Total air temperature in the test section when cloud on 

SH = Specific humidity, g of water vapor/kg of dry air 

Ty = Freezing temperature of bulk water, °C 

Tice = Ice crystal temperature, °C 

Tinelt = Water melt temperature, °C 

T; = Icing surface temperature, °C 

Tail = Wall temperature, °C 

To = Air static temperature, °C 

V = Air velocity, kt 

VWCeain = Vapor water content gained (i.e. = VWC;— VWC;) , g/ m? 

WwW = Fan face mass flow rate, lbm/s 

Po = Collection efficiency at stagnation line, dimensionless 

d = Ice/water mass energy transfer parameter, °C 

Ag Ay = Latent heat of freezing and condensation, cal/g 

7 = Air/vapor energy transfer parameter, °C 

Da = Air density, kg/m* 

Subscripts 

cor = corrected 


ff = fan face 

i = inlet, initial or dry condition 

pl = inthe plenum 

t = target, terminal or wet condition 

S = at the surface 

0 = stagnation value 

1 = PSL-3 contraction duct exit plane, i.e. station | 
a) = local static condition 


I. Introduction 


Mc“ turbofan engine power-loss events have been observed since 1990 on both commuter and large transport 
aircraft near large deep convective weather systems at high altitudes!. The deep convective updraft core can 
bring a large amount of condensed water from the surface to the cloud top. The predicted maximum ice water 
content (IWC) above the freezing level could reach up to about 9 g/m? in theory” at around 30,000 ft. In addition, 
from more recent High Altitude Ice Crystals (HAIC)/High Ice Water Content (HIWC) flight campaign measurement 
data? available on microphysical properties of deep convective clouds it appeared that at cold temperature levels (- 
40°C or below), high IWC areas are primarily composed of small ice crystals (200-400 pm), which yield small 
MMDs around 300 pm. This suggests that high IWC regions at high altitudes occur in and above young convective 
regions where small ice crystal are nucleated in large quantities. 
Based on the instrumented flight-testing data collected in commuter and transport aircraft engine events, Mason 
et al.' has proposed a hypothesis for possible ice accretion due to ice crystal ingestion. The hypothesis recognizes 
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that ice crystals could begin to melt causing mixed-phase conditions within the compression system passage of the 
engine. This mixed-phase water mass would impinge and accrete on unheated surfaces within the engine where the 
local air temperature was significantly higher than freezing before entering the cloud. This ice crystal icing inside 
the engine core flowpath has been shown to create flow blockage leading to loss of thrust and hardware damage due 
to shedding. The hypothesis specifically stated that liquid water is a necessary condition locally for ice accretion to 
continue. Furthermore, it was believed that this ice crystal accretion phenomenon is not unique only to commuter or 
transport aircraft engine type and therefore it was expected that certain common physical mechanisms for ice crystal 
accretion on all jet engines should exist. 

To better understand the fundamental aspects of ice crystal accretion physics and possible influence of high con- 
centrations of ice crystals on engine performance and operability, NASA and National Research Council Canada 
collaborated on a series of fundamental icing physics experimental studies*+>*’ to examine possible physical mecha- 
nisms of ice accretion on surfaces exposed to fully glaciated ice crystal and mixed-phase conditions similar to those 
believed to be encountered during engine power-loss events. 

Key findings from those early studies included: (1) the air wet bulb temperature, T,,,, was identified to be an ef- 
fective parameter to represent the air/water vapor mixture (i.e. humid air) temperature for studying ice-crystal icing 
process, (2) the ice crystal cloud melt ratio, MR, was also identified as an important parameter for the onset and ac- 
tual growth rate of ice-crystal accretion, (3) the ingested amount of IWC was considered a critical parameter directly 
affecting the possible onset and the ice growth rate and (4) there are two distinct types of ice crystal accretion ob- 
served on an unheated icing surface: (i) when the “surface freezing fraction” is positive, it is dominated by the 
freezing of water melt from fully or partially melted ice crystals, the ice structure is formed from accretion with 
strong adhesion to the surface and little water runback, and (ii) when the “surface melting fraction” is positive, it is 
dominated by the further melting of ice crystals, the ice structure is formed from accumulation of un-melted ice 
crystals with relatively weak bonding to the surface and lots of visible surface water runback. 

In 2013 the first NASA ice crystal engine icing test was conducted® on an unmodified ALF502R-5 engine, S/N 
LFO1, which was flown during Honeywell’s 1997 flight test campaign’. This test has shown that the PSL, though 
using the NASA Icing Research Tunnel (IRT) spray bar system technology, was able to successfully replicate sever- 
al known in-flight engine icing events (including the FLT 850 full rollback) and non-rollback events as well. By 
developing such a ground based altitude ice crystal icing simulation facility, detailed aero and thermodynamic pro- 
files of the humid air and ingested ice/water particle size, mass and phase distribution measurement data inside the 
engine core flow path region can be obtained to enable researchers to better understand the underlying physical 
mechanisms in order to further improve and validate the existing in-house engine icing simulation and risk predic- 
tion codes. !*!!!2.!3.4 These engineering modeling tools will be available to the U.S. aviation community to help as- 
sess the potential icing risks and ice protection system requirement as an acceptable Means of Compliance (MOC) 
for icing design and certification purpose. 

In NASA’s current research plan to address the technical challenge of engine icing for the N+2/ N+3 design con- 
cepts there is a milestone to develop an adequate icing scaling method for full engine near sea-level testing which is 
advantageous from a cost saving consideration because altitude facilities are expensive and time consuming to oper- 
ate. The altitude scaling for ice crystal icing is a condition scaling method used to determine the scale test condi- 
tions which would reproduce the same ice shape from the reference condition. Recently a new thermodynamic 
model was developed!> specifically for engine ice crystal icing phenomenon, and the model analysis has been exam- 
ined with aforementioned fundamental ice crystal icing physics experimental data*> and with selected data from the 
2013 full scale ALF 502-R5 LFO1 turbofan engine ice crystal icing test* The preliminary results showed that the 
model is able to capture important qualitative trends of the ice crystal ice formation observed in the engine second 
exit guide vane (EGV2) region under different icing conditions. This model has become the theoretical base of new 
altitude scaling method development!®. The focus of this study is to show some preliminary evaluation results of the 
proposed altitude scaling method from this heavily instrumented ALF502-R5 LF11 engine icing test!’!* conducted 
in PSL in 2015. 


II. Proposed Similarity Parameters for Engine Ice Crystal Icing Scaling Application 


This aforementioned thermodynamic model for engine ice crystal icing (ICI) was primarily based on the scaling 
work originally done by Ruff!? and later extended by Anderson”? and Anderson and Tsao?! using the Messinger’s”” 
steady-state surface energy balance analysis for super-cooled liquid water icing in both Appendix C and super- 
cooled large droplet (SLD) regimes. The expressions for the parameters used will be presented here without much 
discussion. Therefore, readers who are interested in the physical descriptions and detailed derivations of these pa- 
rameters are referred to Ref. 19 and 20 and the references given therein. Several key modifications were made to 
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the surface energy balance equation in order to account for possible interactions between the local air flow within 
the low pressure compressor flow passage and the ingested ice crystals from the engine inlet. 


A. Model Assumptions 

Currently the model only considers the engine ice crystal icing phenomenon for a fully glaciated ice crystal 
cloud which is the most prevalent cloud characteristic in or around a large deep convective weather system at high 
altitudes!. Figure 1 depicts the notional schematics of ice crystal particle path through the ALF502-R5 engine inlet, 
the fan-core rotor and stator (inlet guide vane - IGV), the low pressure compressor (LPC) rotor (supercharger) and 
tandem stators (EGV1, EGV2), the gooseneck duct and the support strut. When the surface energy balance analysis 
is applied at the stagnation region of a LPC component, the model assumes or considers: 

e local air flow is steady relative to accretion process, 


e accretion is quasi-steady, 
e fully glaciated ice particle cloud with single particle size, 
e noLWC supplement (i.e. fully glaciated ice crystal cloud at inlet), 
e when ice particles experience some melting, Tice = Tmen = Ts (i.e. 0 < MR < 1), 
e when MR = 0 (no melting), Tice,i < Tice < Tr, 
e when MR = 1 (fully melted), Ts < Tek < To, 
e for partially-melted ice particle cloud, LWC,; = IWC; * MR, IWC, = [[WC; * (1 — MR)] — VWCegain, 
e onunheated surface (Twan could be fairly hot initially), 
Once ice cloud is on, air starts cooling down, SH LWC,, MVD, {~ local LPC condition) 
increasing and then asymptotic to a steady state IWC,, MMD, 
sueccrrger Melt Ratio (MR) = LWC,/IWG, 
(To, To); Fon Core ‘= a (RH,SH), (Ma, Pols 
(Ma, Po) J <a Support 


(To, Ta) Strut 


(SH, Hy, ; f RE neeeennnse 


e*. Ice Porticle: Poth 7 
WG, Se Lene rntre ar 
Stotor orden Goosenpck << 
MM D; Inet it Stotors Dict 
Inlet /Initial/Dry Target /Terminal/Wet 
(subscript i) (subscript t) 


Melting, Evaporation, Sublimation 
Figure 1. The schematics of ice particle path through an engine fan and LPC region. 
e T,=T; when glaze ice accreted or accumulated on surfaces. 


B. The Surface Energy Balance Equation 

As duly noted in Ref. 4 & 5, there are two distinct types of ice growth mechanism in ice crystal icing. For the 
sake of brevity, only the final equations and key expressions of parameters are presented without much discussion. 
These two types of ice growth mechanism are described as follows: 


(1) Freezing dominated process on surface when cooling > heating (i.e. 0 < no <1) 
conv Bs evap — kinetic =q freeze (1) 


where 
eon = h, (7, =f. -V? /2c 
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pair ), 
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1 i Le is 
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kinetic = Minpt 2 


#nytA, = UIWC, * MR * ff, *V,,) * (1-1, * Mo #A,. 


=m, 


melt imp 


q freeze 


Substituting the expressions for various heat flux terms listed above into Eqn. (1) it can be rewritten in the fol- 


lowing form 
Cc IWC, -VWC,., 
Ny —= D,ws U gain p fe a (2) 
(MR)* A, IWC, b, 


The freezing fraction is formally defined as the ratio of the amount of water that freezes at a given surface loca- 
tion to the total amount of water that impinges at that location. The key terms in this formulation include ¢ and 
6, which have dimensions of temperature and relate to the ice/water mass energy transfer and air/vapor energy 
transfer, and bo, the relative heat factor, which was first introduced by Tribus, et al.” 


Vv 

p =| = ’ 
2c 

Pp, 


Fi 


; ion | eu ie a) 
V ho| T, \t) P 


Ss <2 foe) + me —— Fy fo. ™~ © Vv 
2G h, 1 1" Py us 
0.622\ T, jaa 


(IWC, Z WC ) + B, i (1 - Ht ), - V; . ¢ 


imp,t p,ws p.ws 
by aad =—-_> —, — | rh 
h. h. 
With some algebraic manipulations, Eq. (2) can be further expressed as 
1 * 
Nea Taga G (3) 
by 


in which T\,»,0 denotes the total air wet bulb temperature and is given as 


+6 F oe | (R.H.) 
Ty = T a ea eG: roe) 0 (oo) A ; 


26 x h, 1 le Pry ws ' 
0.622\ T, ie 


and 
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a (oes) IWC, Mn (MR)A IWC, /P, 
0 | ee _vwely ) 


Cp,, IWC, —VWC, ain h, (IWC, —VWC in /P.) 
v2 Y RT. \M? 
oe oh IWE=VWC. |. (Z ay ine es IWC, /‘p, 
=| RA, WC (MR)A, || (IWC,-VWC,,,,)/p, } 


From Eq. (3), it is apparent that both parameters T,,,,0 and MR can affect the accretion process but they are not 
unique. There are many combinations of 7\,»,0 and MR that lead to the same freezing fraction no. It is also noted that 
for a given local thermal profile, i.e. fixed T,,,0 and MR, the freezing fraction no decreases as (IWC;/p,) increases. 
This non-dimensional term (IWC;/pa) suggests that for engine icing problem the air density compressibility effect 
via altitude pressure is also important. For the same amount of IWC; ingested by the engine, the local effective IWC 
value, (IWC/p,), becomes larger with decreasing density due to increasing altitude. This effect was first noticed and 
reported by Veres” ef al. in their recent engine icing analysis work in terms of the ice-water flow rate to air flow 
rate ratio, WAR, from this heavily instrumented LF11 engine icing test data. 

This clearly suggests that the most important parameter affecting the ice crystal accretion process is the freezing 
fraction no which is a function of IWC;/pa), Two and MR. 


(2) Melting dominated process on surface when heating > cooling (i.e. 0 < mo <1) 


kinetic — Tconv ~ evap = FDnelt (4) 
Carrying out the similar substitution procedure as described above, Eq. (4) can be rewritten as 
MR MR a 
6 = Cre fo re Co 
* VWC .,, VWe.,, 
bo 1 —-MR-(——=)) (l—-MR—-(——*")) (5) 
IWC. IWC. 


Similar to the freezing dominated process, Eq. (5) also shows that for melting dominated process on icing surfaces 
both parameters 7\,»,0 and MR can affect the resulting ice accumulation structure through the melting fraction mo but 
they are not unique. In addition, for a given Ty», and MR, the surface ice growth and size (before ice shedding oc- 
curs) will increase as (IWC;j/pa) is increased while the melting fraction mo decreases. 


II. Scaling Method Evaluation from LF11 Engine Icing Test Data 


NASA PSL Ice Crystal Icing Heavily Instrumented Engine Test 
As shown from Eq. (3) and (5), for possible similarity of ice crystal accretion on specific unheated surface loca- 
tion within the engine LPC system flow passage, it is necessary for the scale test to simulate the geometry, the flow- 
field, the ice/water particle trajectories, the total water catch, the heat transfer and, probably, the surface phenomena 
of the desired reference icing condition. However most of the needed information to conduct proper scaling calcula- 
tion for this LF11 engine were not available at the time. Thus an alternative approach to find potential scale test 
conditions was established as follows: 
(1) The FLT 850 case EGV2 icing condition was chosen as the reference condition for this scaling study, 
(2) A number of scale test conditions, with similar engine ICI risk at an altitude of 5,000 ft (5k ft) with best 
matching engine fan-core and LPC operation characteristics (ex: the engine fan face Mach number Me, 
the corrected engine fan face mass flow rate Wor, and the corrected engine fan rotational speed N1 cor) 
from FLT 850, were calculated with NASA in-house codes COMDES” and MELT by matching the 
following three icing parameters of the reference condition at the EGV2 trailing edge and fan face loca- 
tions: 
a. the local static wet bulb temperature T,,o,. 
b. the local melt ratio (WRecvz, re) of a prescribed smaller ice particle size (i.e. 5 um) and 
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c. the engine fan face WAR of the ice crystal cloud (i.e. the non-dimensional ice water content = 
IWC/p,). Ideally, it should match the local WAR if the information or measurements of parti- 
cle size redistribution over the supercharger rotor (see Figure 1) were available. 

(3) A parametric sweep of IWC, MMD, NI and T,) around each proposed scale condition will be performed 
(within the facility operation limit) to ensure sufficient coverage of local reference icing conditions at 
EGV? trailing edge would be simulated to allow reasonable evaluation of such proposed scaling meth- 
od. 

It should be noted that the PSL engine icing test cell continues working on improving its ability to have better ic- 
ing cloud quality in terms of cloud uniformity. The in-house codes are still under development and are very much 
limited by the current state of art capability in instrumentation and measurement to acquire high quality local physi- 
cal variable data inside the engine compression system for ice crystal icing simulation verification and validation. In 
addition, there is no safe and easy close-up visual access to engine LPC exit guide vanes near shroud leading edge 
regions where early ice accretion and shedding were taking place. With those constraints in mind, the images and 
the time-lapsed movies of how the ice structure was formed in the downstream area of EGV2 trailing edge along the 
outer shroud casing surface represent the only available visual characterization of resulting ice accretion. The quali- 
tative comparison of ice formation allows one to visually evaluate how closely this simplified altitude scaling ap- 
proach could simulate the reference ice crystal accretion features. 


A. The Reference and Scale Test Conditions 
The actual run conditions of the reference case (Rdg 93) and two 5k ft scale test cases (Rdg 406 and Rdg 416) 


chosen for comparison are listed as follows in Table 1. 


Table 1. The reference FLT 850 and two 5k ft scale test conditions in 2015 PSL LF11 engine icing test 


PSL ESCORT DATA COMDES/MELT OUTPUT 
Escort | Alt : T. nnvies Twos, | Tsurt, MR RH 
Rdg# (kft) a (ec) To (aim) thks Bae anes EGv2 eev2Te | Wor CRB 
(is) (is) ue (%) 
93 28.9 1.00 0.52 -30.0 1.000 2.3 23 1.00 5.05 -1.32 2.96 0 7.9 1.00 1.06 ae 
406 489 0.95 0.19 -20.8 0.994 5.4 24 1.05 4.51 -1.48 3.68 0 10.6 1.03 1.01 


(close) 


416 486 095 0.19 -26.4 0.971 6.4 18 1.07 5.19 -5.26 1.79 0 10.0 1.05 1.02 pe 


* * * * oe 


Noted in this table that (1) the first eight columns are the PSL designated run number and seven specific meas- 
urement data output from PSL Escort data system, (2) the next eight columns represent calculated quantities ob- 
tained from COMDES and MELT codes and (3) the last column shows the net icing effect (i.e. rollback (RB) or not) 
on this engine and the time taken to call a rollback (CRB). In addition some quantities highlighted with * symbols 
were further normalized by the corresponding reference values, e.g. the engine actual and corrected fan speed (N1’, 
N1 cor), the PSL plenum air total temperature (Tp), the engine fan face Mach number (M“‘) and the corrected fan 
face mass flow rate (W cor). 


B. The Reference and Scale Test Fan-core and LPC Operation Comparison 

In identifying the possible scale test conditions at lower altitudes special attention was given to those with simi- 
lar engine fan-core and LPC operation characteristics to the reference FLT 850 setting. It was noticed from the table 
that although the physical ambient Mach number and fan speed are different, both scale test cases do have similar 
normalized fan face Mach number, the corrected fan speed and the corrected fan face mass flow rate as the reference 
case for this LF11 engine. Furthermore Figure 2 which is a standard engine fan-core/LPC characteristic perfor- 
mance plot also shows that the engine LPC operation setting was fairly close between the reference and those two 5k 
ft scale test cases. 
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COMDES Calculation 


™@ Ref, Rdg93, RB 
@ 5K, Rdg406, NRB* 
A 5K, Rdg416, CRB 


©) 


Pressure Ratio 


Air Flow Rate, corrected 


Figure 2. The LF11 engine fan-core and LPC performance calculation points for the reference FLT 850 and 
those two 5k ft scale test conditions 


C. The Reference and Scale Test ICI Simulation Comparison 

As stated in the beginning of this section, the 5k ft scale test conditions should match the local static wet bulb 
temperature and the local melt ratio of smaller ice particle of 5 um size at the EGV2 trailing edge (TE) location, and 
the engine fan face WAR value of the reference condition to ensure a similar engine ice crystal icing risk based on 
the well calibrated 1-D COMDES code calculation for this engine. However this does not guarantee a similar ice 
accretion growth at the EGV2 TE location and a similar simulation of icing effect on the engine performance since 
the COMDES code currently still does not have a way yet to account for important leading order effects from larger 
ice particle impact, mixed phase ice particle and icing surface interaction, surface runback and surface convective 
heat transfer characteristics. Thus a parametric sweep of [WC, MMD, NI and T,: around each proposed scale condi- 
tion will be performed (within the facility operation limit) to ensure sufficient coverage (i.e. the key is to capture the 
minimum IWC threshold from NRB to RB or vice versa) of local reference icing conditions at the EGV2 trailing 
edge would be simulated. This could provide us some useful information (besides matching YWAR) regarding how 
the scale IWC value at the lower altitude testing should be scaled from the reference IWC value. 


(1) Rdg 93 vs Rdg 406 
Figure 3 depicts a sequence of ice buildup images on the EGV2 TE outer shroud surface area for the reference 


FLT 850 case. The time-lapsed movie suggests the local ice accretion captured by the internal cameras is a melting- 
dominated ice formation process with frequent ice shedding of varying intensity. It should be noted that the refer- 
ence IWC value of 2.3 g/m? has exceeded the minimum IWC threshold needed for a RB event to occur. 


Figure 3. Two ice buildup images on the EGV2 TE outer shroud area captured by the internal cameras for 
the reference FLT 850 case. The time progresses from left to right image, and the flow direction in each 
image is moving from right to left. 


Figure 4 shows another sequence of ice buildup images on the EGV2 TE outer shroud surface area for the 5k ft 
scale test case Rdg 406. The time-lapsed movie suggests the local ice accretion is also a melting-dominated ice 
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formation process with similar reference local static wet bulb temperature and melt ratio. However it became obvi- 
ous during the post data analysis that the scale test [WC value of 5.4 g/m? did not exceed the minimum IWC thresh- 
old needed for a reference like RB event to occur at this altitude. Judging from the EGV2 surface metal thermal 
couple (TC) reading average value over the spray time and the ice shedding frequency and intensity it was estimated 
that if the scale IWC could increase by 20-30% or more, a CRB event with similar reference ice growth features 
would occur. This is an expected tradeoff in using this 1-D COMDES code, though heavily calibrated by the engine 
icing test data, to identify similar engine icing risk conditions at different altitudes for scaling evaluation. It still 
needs a relatively wider range of IWC to ensure the scale condition sweep would capture (or surround) the target 
reference icing condition conservatively. However for this scale test condition the IWC sweep clearly did not go 
large enough. 


Figure 4. Two ice buildup images on the EGV2 TE outer shroud area captured by the internal cameras for 
the 5k ft scale test case Rdg 406. The time progresses from left to right image, and the flow direction in each 
image is moving from right to left. 


After this test point, the LF11 engine icing test team decided to drop the plenum air temperature T,) by 5 °C to 
conduct another similar parametric sweep to see whether a RB event could occur. 


(2) Rdg 416 

Figure 5 displays a couple of ice buildup images in time for this new 5k ft scale test anchor point Rdg 416. 
The time-lapsed movie clearly shows that at this colder static wet bulb temperature and near freezing EGV2 
metal surface temperature the local ice accretion becomes a freezing-dominated ice formation process with firm- 
er ice sheet buildup and less frequent ice shedding. This simulated icing event was occurring at a slower rate of 
thrust reduction even with a much larger scale IWC value of 6.4 g/m? (about 20% more than Rdg 406 case). 
This suggests a higher minimum IWC threshold is required for a RB event to occur in this freezing-dominated 
regime like the Rdg 416.case. 


Figure 5. Two time sequences of ice buildup image on the EGV2 TE downstream shroud area captured by 
the internal cameras for the 5k ft scale test case Rdg 416. The time sequence progresses from left to right, 
and the flow direction is moving from right to left. 
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D. The Reference and Scale Test ICI Engine Response Simulation Comparison 

The reference Rdg 93 case is a full RB simulation of FLT 850 field event in PSL with this heavily instrumented 
LF11 engine. Figures 6-8 show three selected time sequences of ice buildup images along the EGV2 TE outer 
shroud surface area and the corresponding engine responses!” (load, fan and core speeds, booster EGV average met- 
al surface temperatures) during this full RB simulation. Noted that the engine response plots shown here also have 
included the comparison with the LFO1 engine full RB point Rdg 193 obtained from the 2013 inaugural engine ice 
crystal icing test campaign. For readers who are interested in how ice crystal icing can impact this LF11 engine 
performance they are referred to Refs. 17 and 18 for more details. After some careful examination of those figures, 
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Figure 6(b). Engine response at sequence 1 
Rollback Indicators, RD RD 
—— EGVITE RDG1 —— EGV2LE RDG1 ——~ %N2 RDG1 
EGV2TE RDG1 EGV2TE RDG2 ——0S_AVG RDG1 i 
— —%N1RDG1 —load_RDG1 — cloud off_RDG1 
e EE eee S| 
= a 
2 
: 
: : 
z & 
& s 
‘ 2 
Scan (approximately 1hz) ipeeoae 
ROG2 193 
: F 
Figure 7(b). Engine response at sequence 2 
Rollback Indicato RD RD 
—— EGVITE RDG1 —— EGV2LE RDG1 ——~ %N2 RDG1 
EGV2TE RDG1 EGV2TE RDG2 —— OS_AVG RDG1 
— —%N1RDG1 —load_RDG1 — cloud off_RDG1 
a. ES ae le ES S| 
= a 
g 
£ 
: : 
z iy 
& S 
2 z 
‘Scan (approximately 1hz) ern 
ROG2 193 


Figure 8(a). Ice buildup at sequence 3 Figure 8(b). Engine response at sequence 3 


10 
American Institute of Aeronautics and Astronautics 


it becomes apparent that the view captured is not of high value to the scaling method evaluation. Clearly the ice 
buildup on the outer shroud alone could not explain the engine performance degradation seen from those figures. 
There are other critical ice accretion processes hidden away from the camera. Therefore it is imperative to develop a 
better ice shape characterization method to help one understand what critical ice shape features at critical locations 
(could be challenging to have access though) impact the engine performance the most since not all ice accretion 
types pose engine risk. 

As for those two selected 5k ft scale test cases, i.e. Rdg 406 and 416, they were all run and evaluated in a called 
rollback mode. Consequently, there are no detailed comparison of engine response simulation could be made for 
them with the reference Rdg 93 case. A rather simple called rollback (CRB) or no rollback (NRB) and the time to 
CRB were used to assess how well this simplified scaling approach could be used to simulate the reference engine 
response from ICI. The Rdg 406 case was a NRB with a reference like ICI accretion feature. However the scale 
IWC value used was not large enough to have a CRB. As for the Rdg 416 case, it was a CRB, but it was different 
from the reference FLT 850 case. It had a firmer ice sheet buildup and infrequent ice shedding which led to a slower 
rate of thrust reduction and a larger minimum IWC threshold needed for a freezing-dominated RB to occur. As a 
simple comparison, the CRB time for Rdg 416 case is nearly six times longer than the reference Rdg 93 case. 


IV. Conclusion 


A simplified approach for altitude scaling consideration has been evaluated for the 2015 LF11 engine icing test 
campaign conducted in the NASA Glenn Research Center Propulsion Systems Laboratory. This approach first uti- 
lizes a LFOI/LF11 calibrated NASA in-house engine icing risk assessment code COMDES to determine possible 
scale test conditions at 5k ft altitude with similar engine ICI risk by matching the following three icing parameters of 
the reference condition at the EGV2 trailing edge location while maintaining similar fan-core and LPC operation 
characteristics as the reference : 

(1) the local air static wet bulb temperature, 
(2) the local ice crystal melt ratio of a prescribed smaller ice particle of 5 um size, and 
(3) the engine fan face TWAR value of the ice crystal cloud. 

A parametric sweep of IWC, MMD, N1 and T,) around each 5k ft scale condition was performed to better under- 
stand both the ice accretion and the engine icing effect characteristics. This would help evaluate how similarly the 
reference ICI and engine icing effect features being simulated. 

The instantaneous images and the time-lapsed movies of ice buildup in the EGV2 TE region could only provide 
limited ICI accretion characteristic information. A better ice shape characterization method with respect to the ice 
shape effects on engine is clearly needed. Those ice buildup images on the engine EGV2 outer shroud area are used 
with other relevant air/ice crystal thermodynamic and engine performance data to qualitatively and visually assess 
the 5k ft engine ICI simulation results. It is found from the preliminary evaluation that altitude scaling is possible 
for PSL or a sea level engine test facility to simulate the important ice crystal accretion features and the associated 
ICI effects on engine performance. But special considerations and unique test methods are needed in order to ad- 
dress the associated facility operation limitations, e.g. altitude vs sea level or direct connect vs blow through test 
facility, for running different type of engine icing test configurations, for example a full engine vs a motor driven 
rig. 
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